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Abstract 

We examine the interplay between size quantization and Landau quantization in the De Haas- Van 
Alphen oscillations of clean, metallic nanowires in a longitudinal magnetic field for 'hard' boundary 
conditions, i.e. those of an infinite round well, as opposed to the 'soft' parabolically confined 
boundary conditions previously treated in Alexandrov and Kabanov (Phys. Rev. Lett. 95, 076601 
(2005) (AK)). We find that there exist two fundamental frequencies as opposed to the one found in 
bulk systems and the three frequencies found by AK with soft boundary counditions. In addition, 
we find that the additional 'magic resonances' of AK may be also observed in the infinite well case, 
though they are now damped. We also compare the numerically generated energy spectrum of the 
infinite well potential with that of our analytic approximation, and compare calculations of the 
oscillatory portions of the thermodynamic quantities for both models. 

PACS numbers: 75.75.+a, 73.63.Nm, 73.63.b 
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I. INTRODUCTION 



Since their discovery in 1930 [lj|, and subsequent explanation by Onsager 2] and Lifshitz 
and Kosevitch jj] in terms of the interplay between the Landau quantization and the Fermi 
surface, the de Haas- Van Alphen (DHVA) oscillations of the thermodynamic potential, mag- 
netization and related quantities has acquired a long history of use as a method of probing 
the structure of the Fermi surface of metals at low temperatures and large magnetic fields. 
(A summary of the first thirty-eight years of study in this field may be found in Gold J4]. 



More contemporary reviews of the area may be found in Shoenberg 5[, Singleton [6j and 
Kartsovnik 7], with the latter two focusing on applications to quasi-2D organic compounds.) 

The question of boundary conditions and how confinement affects the behavior of elec- 
trons in longitudinal magnetic fields (and hence the DHVA effect) has long been a question 
of interest to the condensed matter community (see, for example, Fock jsj], Landau 9], Dar- 
win [h]] and also the discussion in Dingle's classic paper on the effects of boundary 
conditions on the density of states in large systems). This issue is of obvious importance 
in the study of nanowires, where the system of interest is confined within a narrow, cylin- 
drical potential, which in this case is difficult to treat analytically owing to the difficulties 
in approximating the eigenvalues of the Schrodinger equation given by the zeroes of the 
confluent hypergeometric function. However, in the weak field limit, the eigenvalues are well 
approximated by the zeroes of the Bessel function, and so Dingle 12| was able to calculate 
the effect of such a confinement in a weak magnetic field, through treating the latter as a 
small perturbation to the energy levels and showing that at low fields, one should expect 
to observe oscillations that are primarily due to the passage of the size quantized energy 
levels above the Fermi energy. Bogachek and Gogazde 13] calculated the contribution to 



the thermodynamic potential of those electrons near to the surface of the wire that are on 
a 'grazing' trajectory with respect to the surface, and argued that this is in fact a more 
dominant contribution to the thermodynamic potential than the Dingle-type oscillations, 
as the latter are suppressed by specular reflection under normal circumstances. The con- 
tribution of these surface levels corresponds to that induced by the Aharonov-Bohm effect. 
This effect was observed, along with others ascribed to Dingle-type oscillations by Brandt 



et al. 



14[ in a study of the magnetoresistance of bismuth wires of micrometer thickness. 



Further confirmations of the presence of the low field Aharonov-Bohm effect before the on- 
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set of Shubnikov-De Haas oscillations in the magnetoresistance comes from the study of 
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progressively narrower Bi nanowires [15|, as well as in carbon nanotubes 

It is likely that at higher fields for wires of very narrow radius, we might expect addi- 
tional effects of the interplay between size quantization and Landau quantization to become 
apparent. There is a need for theoretical treatments of this region, and in what follows we 
make one of the first attempts at doing so outside of an idealized 1-D limit, such as that 



used in Sineavsky et al. 



Alexandrov and Kabanov (AK)[19] have recently examined the DHVA oscillations present 
in a system of electrons confined within a parabolic confining potential (the Fock-Darwin or 
confinement model 0, [h]]). They have found that there exist three characteristic frequencies 
of DHVA oscillation, and also that the thermodynamic potential displays 'magic resonances' 
at certain values of the field where the amplitude is greatly enhanced. In a brief note at the 
end of this Introduction, we demonstrate that despite being due to a divergent cotangent 
function in our expression for the oscillatory thermodynamic potential Q, the behavior of 
the thermodynamic potential is still analytically well defined - this is important if the AK 
result is to be properly understood. 

Our paper focuses on a system of free electrons confined within a cylindrical infinite 
well potential and subject to a longitudinal magnetic field. Using a well known asymptotic 
description of the wavefunction of such a system, we derive a non-linear energy spectrum and 
confirm it by deriving a similar result through a semiclassical approximation. By linearizing 
these spectra near the Fermi energy, we find that we acquire expressions for the energy 
similar to that of the Fock-Darwin model, and so we may motivate an improvement to 
the asymptotic and semi-classical models that brings them closer to the energy spectrum at 
larger values of the field than was previously possible. The oscillations in the thermodynamic 
potential are calculated, and we find not only that the additional characteristic frequencies of 
the system are modified, but that the position of the magic resonances has altered and that 
they are somewhat damped. Numerical plots of the thermodynamic quantities calculated 
are given for various values of the field. 

Finally, we examine how well the modified asymptotic approximates the exact spectrum 
of the problem. We find that it is a reasonable approximation for low values of the magnetic 
quantum number m at large u c = eB/m*c; however, it is inaccurate for large values of m, 
reflecting orbits that are closer to the edge of the wire. This is likely because the asymptotic 
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that we are using is not valid in that limit (similarly, the semiclassical approximation breaks 
down here since it is too close to a turning point of the classical electron trajectories). 
However, the actual boundary conditions of the system are most likely intermediate between 
those of the infinite cylindrical well and those of the Fock-Darwin confinement model, and 
so it is perhaps reasonable to suspect that the breakdown of the approximation, which is 
likely due to the 'hardness' of the boundary conditions in the infinite well case, should not 
unduly restrict the application of the nonlinear spectrum to real systems. 

In what follows, we assume that the system of interest is a long, clean metallic nanowire 
in a longitudinal magnetic field B that is parallel to the direction of the wire (taken to 
be the z axis of a cylindrical coordinate system). We take the electron mean free path, 
I = vpT to be comparable to or larger than the radius of the wire R, but not the length 
L. We also assume that the electron wavelength near the Fermi level is very small in these 
metallic nanowires, such that L ^> I < R ^> 2irh/(m*v F ); here and previously, vp is the 
Fermi velocity and m* is the band mass in the bulk metal. 

From henceforth, we also take H — c — m* — L — 1. 



A. On the divergence of the 'magic resonances' in the AK result 

The AK [19] result describes the DHVA oscillations of a three-dimensional electron ex- 



posed to a magnetic field within a two-dimensional parabolic confining potential [8|, Il0 |. 
This system is an isotropic planar harmonic oscillator whose frequency in the absence of a 
magnetic field is taken to be ujq. Due to the effects of confinement, the DHVA oscillations 
of the system aquire three characteristic frquencies of oscillation: 2u, uo + and u~ , where 
oj 2 = Uq + us 2 / 4 and w ± = uj ± to c /2. The amplitudes are given by 

T{2xf/ 2 /irry\ 
r(X ' V) ~ 27rr3/2 s inh[7r 2 Tr/(x)] COt \ 2x ) ' [) 

where x = uj and y = cj ± for the the amplitudes of the oscillations with frequency 2u, and 
x = uj^/2 and y = 2uj for those with frequency u ± . As has been previously noted, these 
amplitudes (0Q) contain a cotangent function which gives rise to 'magic resonances' for certain 
values of the field where the condition 2uj/{uj ± ) = (q + 2)/r is satisfied. If this is the case 
then the cotangent becomes infinite provided that q is an integer. 

However, despite what one might naively assume, this apparent divergence is in fact 
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analytically well denned. Let us examine the following harmonics: r = 1, q = 1 and 
r = 3, q = 9. From our condition for the divergence of the cotangent, we have to = |u; c and 
hence to~ = Let us set u = 1, and then set uj~ = |(1 — 5), where 5 — > +0. 

For our r = 1 harmonic, we shall examine the Ai(^,2u) amplitude, and for our r = 3 
harmonic, we shall examine the A^(uj,uj~) amplitude. We acquire: 

Al (lP 2U ) = Al (^ 1 - 5 )' 2 

T(2/3)V2 (2) 



27rsinh[37r 2 T] V 1 " 5 



and 



4,(0;, or) = ^A^l-*)) 
T(2/3) 1 / 2 



. cot(7r(l-5)). (3) 

07T sinh|37ri J 

Expanding out the cotangent from (j2J), we find cot (37r/(l — 5)) ~ 1/(3tt8) while doing the 
same for (j3J) gives cot (it(l — 6)) ~ — l(7r<5). Substituting these values back into (j2J) and (j3J) 
respectively, we find that: 

We may now write, where Q res is the oscilliatory component of the AK result near the 
magic resonance: 

n res = A 3 (w, uj sin y — L - - + A x — , 2u sm i 

\ uj 2uj 4/ \ 2 J \ uj- uj~ 

A . ( 3tt / 2 M \ rr\ A . ( 3/itt 3tt / 2 M 

Defining C = 3/i7r — |vt(c<j + + 1) and expanding out the terms in the denominators, we arrive 
at: 

Q res = AC cos (C - |) • (5) 

Note that the singular terms are no longer present, so that the thermodynamic potential in 
the region of the resonance is now finite. However, it has also been enhanced by a factor of 
C, which is of the order of the Fermi energy /i ^> 1, and so remains rather large. We would 
expect the same analysis to hold true for the other magic resonances. 



II. CYLINDRICAL CONFINEMENT AND THE DHVA EFFECT 



The Schrodinger equation in polar coordinates for an electron moving in a longitudinal 
magnetic field B parallel to the z axis is: 



1 d_ 

r dr 



dr 



d 2 ib 1 d 2 ib 
H H - 

dz 2 r 2 deb 2 



itu c dip uj 2 r , _ , 



(6) 



Following for example (though note the difference in convention regarding e) Landau and 
Lifshitz [20] we seek a solution in the form 

f 



R(r)e tm(p e ip * 



where R(r), after redefinition in terms of £ = (u c /2)r 2 is the radial function: 

L 



|m| 

IT 



\m\ + l,£), 



(7) 



(8) 



H 
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I). 



m| + 1, £) is a confluent hypergeometric 



with/? = tu- 1 [E-(p 2 z /2)]+m/2. Mal- 
function. The zeroes of this are the eigenvalues of the Schrodinger equation for cylindrically 
confined electrons. If we write (3 — ^ — | as —a, followin g th e notation of Abramowitz and 
Stegun 2l|, the energy levels of the particle are given by 20] 



k 2 

E = h uo c 

2 



-a + 



m + 1 



(9) 



If the wave-function is finite everywhere, then —a is an integer, and the eigenfunctions of 
the equation are the Laguerre polynomials. However, if this is not the case, and it becomes 
zero at some finite radius R, then —a is a positive, real number. In order to proceed 
analytically in these cases, particularly where R is small, we must therefore use various 
approximation or numerical methods so as to obtain the energy spectrum. 



A. Asymptotic approximation 



For large enough values of a, where a is negative, the confluent hypergeometric function 



M(a, b, £) may be approximated by the following asymptotic form 
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M(a, 6, 0=r(6)e«/ 2 



(l-6)/2 



7T 



"3 cos ^a/26£ - 4< - ifor + 



(10) 



where (as before) we take a = —(uj c 1 [E — \ + y — ^ — |), b = \m\ + 1 and £ = u c R 2 /2 = 
(-7T /2ui s ) 2 imjj c1 with i? being the radius of the nanowire and u s = ixv F /2R. In so choosing 



the value of £, we are imposing 'hard' boundary conditions by stating that the wavefunction 
must be zero at R. 

Our approximation has zeroes when the following is true: 

V2&e - 4< - \bn + ±tt = ^±lvr, (11) 

and so, after some algebra, we obtain the following approximate expression for the eigenval- 
ues, valid wherever the asymptotic ( fTOl) is valid: 

2 / II o\ 2 

lo I \m\ 6 \ uj c m 



A similar energy spectrum may be derived from the Bohr-Sommerfeld quantization con- 
ditions: 

2 / ~ , 1 \ 2 

oj. I m + 1 \ u; c m 



E = 7f(" + — J "i-" (13) 

where m = -^/m 2 — 1/4. 

Up to a phase factor this is very similar to (fl2|) . which is derived from the asymptotic 
form of the confluent hypergeometric function. In what follows we shall tend to use the 
spectrum (fl2i) unless otherwise stated, since there is a possibility that the Bohr-Sommerfeld 
conditions do not work as well as they might, as has been noted in the case of a nanowire 
in a transverse electric field (as observed in 221], for example). 



B. Connecting the approximations with the parabolic confinement model 



We now turn to address what, if anything, (Tl2l) and ( fl3l) have to do with the confinement 
potential spectrum ( jg, 10, 3] - here we take u to be the frequency of the confining potential 
used in those sources). For simplicity, we shall set u c = (i.e. work in the limit of zero 
magnetic field). Defining £ = E — fi, and working where E w fi, we may write for our 
non- linear spectra (11211131) : 

(14) 



£ = ^A 2 -fim 2u s A - 2fi 



where 



A 



i \ m \ i 3 

n + + I 



asymptotic 



(15) 



[n + semiclassical 
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When lo c = 0, E « fi, the spectrum of the parabolic confinement model becomes: 

/ 1777,1 + 1 \ 

£ = E-ii = 2u [n+ l —±-—) -ii. (16) 



So, in this limit, the confinement model, if we take ujq = u s , and the non-linear approx- 
imations are good matches up to an unimportant phase and a factor of 2 in front of the \i 
term. Since the latter only shifts the value of the energy by a constant, it seems that the 
results of are a good approximation of the behavior of the non-linear approximations 
( fl2lfT3l) . very near the Fermi energy, and at zero field. 

This motivates a replacement of u s in (Tl4|) . and hence in ([121 and (fl3|) . with uo so as to 
better approximate the behavior of the system at larger values of uj c with these non-linear 
spectra. That this is not a bad approximation (at least for small and intermediate values of 
m) is suggested by the numerical analysis of Section HVl 

Having made this replacement and restoring the momentum component in the z direction, 
we may write the spectra of the non-linear model ( 1121131) of electrons confined in a nanowire 
as: 

^ k 2 to 2 , 9 uj c m . . 

E= 2 + 7*-—- (17) 

It is also instructive to compare this with the spectrum used by Bogachek and Gogadze 
[13I ] and Dingle jl2| . Perturbation theory for small values of the field gives us the following 
expression for the energy: 

f k 2 

2m* R 2 2m*' 1 ' 



where in this case we restore m* for clarity, and 



1 



2(m 



2 



I 2 

1 nm 



(19) 



with 7„ m being the n-th zero of the m-th order Bessel function, and <fi = $/$o? where 
$ = irR 2 B and $0 = hc/e (in ordinary units). Defining a frequency v a = vp/R our overall 
expression for the energy becomes: 



k 2 v 2 2 uj c <fiu c 



E = — + - — + 



1 + 



2(m 



2 



7 2 
Inm 



(20) 



2 Afi inm 2 12 

Discarding the last term, and noting that the zeros of the Bessel function in our region 
of interest (as opposed to that of Bogacheck and Gogazde 13] and Dingle [12]) are given 
by n(n + ^ + ~) we arrive at our asymptotic expression for the energy at weak field. It 
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is possible to recover the Aharanov-Bohm oscillations predicted by Bogacheck and Gogazde 
from the results of our calculations in the next section by examining the limit of small u c . 
However, the oscillations they obtain have the form cos(27rr0), whereas ours have the form 
cos(4r0/7r). It seems probable that for the most part this discrepancy due to our use of a 
different asymptotic approximation of the Bessel function zeroes - they use an approximation 
valid for large m, and we use one valid for large n. 



C. The De Haas Van Alphen effect for the non-linear model 

We now calculate the oscillatory portion of the thermodynamic potential for the asymp- 
totic version of (fTTl) - henceforth referred to as the non-linear model - neglecting the phase 
term in the quadratic portion since it is much smaller than n + ^p, so that we have: 

j 2 2 

E nmk = y~V + 7 (n + ^ ■ (21) 

Applying twice the Poisson's formula to the expression for the thermodynamic potential 

n = -T^ln[l + e ( ^ )/T ], (22) 

a 

and replacing negative m with — (m + 1) one obtains 

oo 

n = E F& + n^]> (23) 

r,r'=— oo 

where 

95/27-1 5/2 /-oo poo poo 

n±, = - - -4 / dk dx d 2/e ^+w)ln[l + e ^ 1 -(^) 2 - fc2± «/ T ], (24) 
^ Jo Jo Jo 

p = 2TCfir/u), q = 4irfir'/u}, and j3 = uj c /uj. We are interested in the part of Q oscillating 

with B, which arises from terms in Eq.(4) with nonzero r,r'. Replacing (x + y) for z and 

introducing the polar-spherical coordinates yield 

95/2T-1 5/2 poo poo 

= / dye l{q - p)y d<f> dppe ippcos ' t> \n[l + e^ 1 - f,2±M/T }. (25) 

nuJ Jo Jo Jy/coa<f> 

Replacing y for y = (p 2 — e)/ (±/3) we obtain 

95/271 5/2 poo p-w/2 pl{e,4>) 

™ 2 (3 J Jo J^i 

(26) 



where j(e, 0) — [e + f3 2 cos 2 (0)/4] 1//2 + /3cos(0)/2. f2~, is obtained by replacing (5 for — /3 
in this expression. Here we neglect a contribution from regions where the total energy e is 
negative. This contribution is an artefact of the approximation (12~T1) . and, in any case, it is 
exponentially small at temperatures T <C p in the oscillating part of Q. 

Neglecting terms of the order of u/p <C 1 the integral over p is approximated as 

j(e,<t>) \ f- ipyjl cos (j,+i(q-p)e /P f A\ p ipy(e,<j>) cos cj>+i{q-ph 2 (t,<l>) / P 

^ p(3 cos <fi + 2(q - p)e p(3 cos <fi + 2(q - p)^{e, 0) 

(27) ■ 

At low temperatures, T <C /i, the main contribution to the oscillating part of the thermo- 
dynamic potential, Q, comes from energies near the Fermi surface, |1 — e| <C 1. Moreover at 
large p ^> 1 only small angles <C 1 contribute to the integral, which allows us to extend 
the integration over up to infinity. We can also replace cos(0) by 1 — 2 /2, expand 7(e, 0) 
as 

^ K -<-wm~wm' (28) 



in the exponents (here 7 = ^/l + /3 2 /4 + /3/2), and take cos0 = 1, 7(e, 0) = 7 in the 
pre-exponential terms. Integrating over and [ 1 — e | with the use of J °° o?0exp(ia0 2 ) = 
{TijA\a\) 1 / 2 exp[z7ra/(4|a|)] and J dxexp(iax) [1 + exp(a;)]~ 1 = —277/ sinh(7ra) yields 

e ip-iirp/A\p\ 



~ , 87T 1 /2 T/ i 5 /2 



^ 2 



plpl 1 / 2 sinh(7rpT/2/i)[p/3 + 2(g — p)] 

(ry _ ^^2) 3 / 2 e i97_i7r9 / 4 l 9 l 



(29) 



g|g|V2 s inh[7rgT/2/i(7 - (3/2)}[p(3 + 2(g - p)j] _ 

This approximation fails at the 'magic' resonances, where pft + 2(q — p) = or pj3 + 2(q — 
p) r y = 0. In those cases the integral (1271) . is calculated as 

/ PS (7T/4|p| ) 1 / 2 e i F'V / e cos p)e//3— -ivrp/4|p| ^ j-gg^ 

or 

/ ^ 7(7r(7 — /3/2)/4|g|) 1 / 2 e ip7 ^'^ cos ^ + ^ g_p)72 ^'^/ /3 ~ i7r ^ 4 ' 9 ', (31) 

respectively. Combining equation (1281 and equations (130113 X I) we can replace the resonant 
denominators in equation (|2"T|) as 

1 1 



p(3 + 2(q-p) p(3 + 2(q-p) + i/3(4|p|/7r) 1 / 2 exp(i7rp/4|p|) ' 



(32) 
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and 

1 1 

pP + 2{q - p) 7 ^ pP + 2{q - V )l ~ i&W\Kl ~ P/ 2 )Y /2 exp(mg/4|g|) ' 
Substituting this expression into ( l23l) and performing partial summation yields 



where 



and 



(33) 



VI = 5^^^ sin(27rr/i/cj - 7r/4) + Bf cos(2irrfJ,/u - it /A) 

r=l ± 

+ Cf sin[47rr / u/d; =F - tt/4] + Df cos[47rr/i/d; =F - vr/4] (34) 



A - = 2^SI&M Scot|TrwT/w * 3(1 " 'N™'"^ (35) 

B ? = 2^1t(^) aCOt ^ T ^ ± 3(1 - 0(~//.W4], (36) 



^ = I 3/2 T< & tT/}- + x ..,,K«»t[^/(^4ai-)±^l- Vfr/rt^/Z^-fui-n, 
ginh[47r / i r/{oo + + u )\ 

(37) 



Anr 3 / 2 smh[47r 2 Tr/(a; + + 



(38) 

Here uj ± = lu±lu c /2, = + (3 2 /4±(3/2), and the summation formula, ^ r (-2 — r) _1 = 

7rcot(7rz) has been applied. The +/— of Af cancel with each other, as do those of Bf - 
this allows us to simplify our formula and write: 



oo 



Q = 2^ 2^ Gf sin[4vrr/i/w T - vr/4] + Df cos[4vrr/i/cD T - vr/4]. (39) 

r=l ± 

We can see that, unlike the AK result, we have two characteristic fequencies u /2 which 
differ from the equivelent frequencies in that result. In addition, the magic resonances occur 
at different ratios of u c / u s and are damped by the additional terms in the cotangent function. 
As is usual in these calculations, we may restore the effects of spin-splitting by multiplying 
each term by cos(riigfj, B /e). It should be noted that the effects of spin orbital coupling and 
other complications may further complicate this expression in real materials. 
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III. NUMERICAL RESULTS 



Here we present some numerical calculations of the DHVA oscillations and of their Fourier 
transforms so that the AK and nonlinear results might be compared. We have set \x = 2000co> s 
and u s = 1, measuring u c and all other quantities in units of uj s . Fourier transfromations 
were performed using the NAG DFT routine C06FAF and the Fourier integral calculation 
techniques described in {23]. 

Calculations at T = for the parabolic confinement model were performed through 
taking the limit T —>■ in (122]) and then integrating over k. We thus obtained the following 
expression for the thermodynamic potential, 

n=-^Y j ^-E nm )l, (40) 

n, m 

where E nm is energy for a given n, m, The non-oscillatory portion was then subtracted out, 
leaving only the oscillatory portion of the function behind. 

This could not be done in the case of the non-linear model, however, due to the issue 
remarked on in Section III CI regarding the unphysical negative energy levels. It appears 
that as yet there is no straightforward way of calculating the thermodynamic quantities 
in this fashion, as it does not appear as easy to remove the contributions from negative 
energy values in this method as it was in the aforementioned analytic ('Poisson summation') 
treatment of Section III CI 

For T > 0, the calculations were done using the formulae of [19| for the AK results, and 
fl39|) - fl38|) in this paper were used for the non-linear results. 

Figure [1] displays the behavior of the thermodynamic potential with respect to changes 
in (cg> c ) -1 , alongside their Fourier transforms. These last clearly show the three peaks cor- 
responding to each of the characteristic frequencies given in [l9|] (here, Fq is equivalent to 
the F of that reference). In the case of the T = plots it seems that some residual non- 
oscillatory components of the untransformed functions remain at low frequencies, and these 
obscure the low frequency content of the Fourier transform. 

In Figure [2] we present the first magic resonance in the thermodynamic potential of the 
non-linear model together with that of the parabolic confinement model for the purposes 
of comparison. We can see that the resonance in the non-linear case has not only changed 
position, but has also been damped by the additional terms in the cotangent functions of 
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FIG. 1: Graphs of Q in the parabolic confinement model opposite their Fourier transforms. Note 
that the remnant of the non-periodic component of the T = data is most noticeable in the 
topmost graph, but also that the F+ resonance is most obviously visible in its Fourier transform. 

the Fourier coefficients. 

Figure [3] shows the thermodynamic potential and its Fourier transform in the same ranges 
as previously used for the parabolic confinement model calculations. Following , we may 
estimate the characteristic frequencies as being: 



2F|\/2± y/T 



7 



(41) 

ttR 2 is the cross- 



v /r+Y(v / r+7 7 ±(i/v / 2)) 2 ' 

where F = S F /{27ce), 7 = 4u 2 s /u 2 c = 4vr 2 S F / \e 2 S B 2 ) , 7' = 7 /2, S 
sectional of the wire, and Sf 
surface. (We remind the reader that in order to interpret the frequencies given in the figures 
using these formulae, we should take e = m* = u s = 1). 



7r(m*) 2 v F is the cross-sectional area of the Fermi- 
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FIG. 2: Graph comparing the first magic resonances of the parabolic confinement model with the 
non-linear model. The temperature here is set at T = 0.025. 

We can see the two characteristic frequencies clearly in the Fourier transform, and that the 
oscillatory behavior of the function is considerably different from that seen in the parabolic 
confinement case. 



IV. HOW GOOD ARE THE APPROXIMATIONS? 



In order to understand how our non-linear model differs from the exact behavior of the 
system, we must compare their spectra with that generated from the exact zeroes of the 
confluent hypergeometric function M(a,b, £). In order to obtain the latter, we numerically 
generated M(ct, 6, £) with the CONHYP function 24], |25|, and located the roots of the 
function as the parameter —a was increased from zero using the Van Wijngaarden-Dekker- 
Brent algorithm 26]. 

Unfortunately, the processing ability of the computer limits the values of the variable <; for 
which results can be obtained; this means that it is difficult to obtain zeroes for large values 
of both fi/u) s and uo c /uj Si since this entails a large value of £ that will cause the series defining 
M(a, b, £) to be very slowly converging. However, one can gain an idea of the behavior of 
the exact system for less extreme values of the variable, and it is to such an example that 
we now turn. 
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FIG. 3: Graphs of 0, in the non-linear model opposite their Fourier transforms. Note that the F+ 
and F— peaks in the lower Fourier transform are so close as to be indistiguishable; this may be 
why the second harmonic is larger than expected. 

We set /i = 100u s , u s = 1 and examine values of < oj c /uj s < 1. For each of these, we 
determine the largest value of n in the non-linear approximation for which E < \x (assuming 
k is zero) and calculate the relative error A = (E — E exact )/E, where E exact is the exact 
energy calculated from the n-th zero of M(a, b, £). This enables us to gain some idea of the 
size of the error at the Fermi surface described by the non-linear model. 

To obtain the exact energy in the zero-field limit, we note that the wavefunction in that 
case takes the form of a Bessel function, and so the exact zeroes may be obtained from the 
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FIG. 4: The relative error A along the Fermi surface defined by the non-linear model, as a function 
of m and uj c . 

function J\ m \{2yJE exact w), where w = (^V 

We plot A in Figure HI It is apparent that as u c increases, the shape of the approximate 
Fermi surface alters such that the maximum possible positive value of m increases, and the 
maximum negative value of m decreases. This coincides with a decrease in the maximum 
value of | A | for negative m, and an increase for positive m - in fact the latter error is so large 
for large m at u c = 1 that it can scarcely be said that it accurately models the exact values 
at all in that region. In addition, it is worth noting that the non-linear model underestimates 
the deformation of the Fermi surface when u c is large (for example, see Figure Why is 
this? 

In the case of the asymptotic derivation of the model, the asymptotic from which we 
obtained it makes two assumptions: that —a is greater than £, so that M(a, b, £) resembles 
a Bessel function, and that —a is greater than b, so that the Bessel function resembles a 
cosine function. Obviously, this condition is violated in the zero-field limit at large m, and 
so we would expect it to be violated for large m at finite values of the field, even though it 
appears that the replacement of u s with to in Section III Bl counteracts the violation of the 
first condition for large £ if b is small. This is what occurs in the case of positive m. As 
for the the portions of the spectrum with negative m, since these have larger energies than 
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FIG. 5: Graph showing the quantum number n against m at the Fermi surface for fx = lOOus and 
uj c /uj s = 1.0. 

their positive counterparts and since this disparity increases as we increase u c , one should 
not be surprised to find that, as, u c increases, the Fermi energy is exceeded at smaller and 
smaller values of negative m, which are in turn more closely matched by the asymptotic. 
And so in this case the maximum error decreases. 

Semiclassically, one should note that large values of m describe electron orbits that ap- 
proach the edge of the wire quite closely. Since that is a turning point due to the imposition 
of the 'hard' boundary conditions, we cannot expect the semiclassical approximation to hold 
in that region. The analysis of the difference between the behavior at positive and negative 
m given above also holds here. 

Perhaps we should consider the boundary conditions of our problem. The parabolic 
confinement model was originally derived with 'soft' boundary conditions - a parabolic 
confining potential - and as we have shown, this approximates the non-linear model, derived 
with 'hard' boundary conditions near the Fermi surface. However, this latter model does 
not capture the large m behavior of the system very well, in part, one suspects, because of 
the 'hardness' of the boundary conditions. It should be noted then, that neither of these 
extremes is likely to be physical - the boundary conditions for a real nanowire will most 
probably be intermediate between them, and so it is possible that the experimental behavior 
of the system might also be intermediate between the behaviors we have discussed. This, of 
course, is a matter that requires empirical determination. 
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V. CONCLUSIONS 



We have calculated the oscillations in the thermodynamic potential due to the De Haas- 
Van Alphen effect in a clean, metallic nanowire with a simple Fermi surface for a nonlinear 
energy spectrum derived from an asymptotic approximation to the exact solution of the 
Schrddinger equation, and compared it with the AK results derived using the parabolic 
confinement model which approximates it near to the Fermi energy. In both models, one 
can observe 'magic resonances' where the amplitude of the oscillations in the thermodynamic 
potential is enhanced at particular ratios of u s to u c , although their locations differ in each 
model, and are somewhat damped in the non-linear case. However, the nonlinear result 
lacks one of the frequencies of oscillation observed in the AK result. In addition, it seems 
that the infinite well boundary conditions cause some problems at large values of m due 
to the semiclassical nature of our approximations; however, this is unlikely to be fatal to 
our predictions of magic resonances since the real system likely has boundary conditions 
intermediate between the two cases which we have considered. 

It seems clear that the interplay between the size quantization and the magnetic quan- 
tization is non-trivial in the extreme, and exhibits a certain degree of model dependence. 
Further experimental and theoretical work is needed both to clarify remaining ambiguities, 
and to apply this general theory to the more complex band structures of realistic nanowires. 
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